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We investigate perturbative aspects of gravity with a general F(R) Lagrangian, as well as non- 
perturbative dilatonic solutions. For the first part, we are interested in stability and the definition 
of asymptotic charges. The main result of this study is that, while generic F(R) theories are stable 
under metric perturbations, they are expected to show instabilities against curvature perturbations 
when the Lagrangian includes 1/R terms. For the second part, one is interested on exact solutions, 
and we explicitly construct kink-like solutions of the Liouville type for the dilaton field for F(R) 
■ having the form R + fR n , in two and in four dimensions. 
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fir I. INTRODUCTION 

< 

. In this work we investigate perturbative and nonperturbative effects in gravity, motivated by the cosmological 
constant problem As one knows, the cosmological constant engenders an important but very hard problem to 
' solve due to both cosmological and phenomenological reasons. In the cosmic evolution of the Universe, particle 
^ ■ physics introduces the main mechanisms to govern the several phases of the Universe, and cosmology drives the way 
the cosmological constant is accounted for in each phase. However, the presence of spontaneous symmetry breaking, 
which is the basic mechanism underlying unification of the fundamental interactions in nature, inevitably injects 
vacuum energy density even if no such term is present in the theory. For this reason, one realizes that not only the 
physics at the Planck scale, but also the low energy infrared dynamics of gravity itself is essential to the resolution 
of the cosmological constant problem. The reasoning establishes direct connection between nonperturbative effects in 
gravity and the cosmological constant problem. 

Nowadays, the cosmological constant problem must include the important discovery that the Universe is currently 
Qh, evolving in an accelerating phase Q|. There are several distinct possibilities of taking acceleration into account: from 
one side, we can keep standard geometry, incorporating modifications on the matter contents of the model. Interesting 
investigations which deal with this possibility include, for instance, the cosmological constant Ij, dynamical scalar 
fields Chaplygin fluid 0], and phantom dynamics Q. Another possibility deals with modifications of geometry, 
changing General Relativity; here an important class of theories includes models which depend on the Ricci scalar, 
usualy named F(R) theories [6]. In the present work, we will deal with the latter, that is, we will investigate extended 
gravity, changing R — > F(R), for F(R) — R + jR n , including the cases n — 2 and n = — 1. This modification brings 
interesting novelties, as we will show below. The motivation goes beyond the recent discovery of the accelerated 
expansion of the Universe Q, and it is very natural to expect that, due to the fact that the curvature of the 
present Universe is very small, the nonperturbative effects should be characteristic for some typical forms of F(R). 
In particular, the interest in n — 2 also appear from semi-classical investigations, and the case of F(R) = exp(— aR) 
suggests the presence of instanton effects. 

The investigation of nonperturbative solutions is of direct interest to a diversity of applications, including black 
hole, integrability, noncommutative geometry and other problems (9I. [Toj| . The nonperturbative issues are in general 
hard to solve, mainly in high energy physics in curved space-time. For this reason, we are planning to approach the 
problem with the use of the nonperturbative procedure referred to as the first-order formalism, in which one solves 
equations of motion with solutions of first-order differential equations. Some of us have recently done some progress 
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in extending the formalism to curved space-time [ll| , and we now consider this possibility in connection with dilaton 
gravity. 

In the present study we organize this Letter as follows. In the next Sec. [TT] we derive the equations of motion 
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for F(R) gravity theories and the conserved charges, and then we study stability of solutions in these theories. In 
Sec. IIIII we turn attention to gravity in the conformal sector, and there we describe the F(R) model restricted to the 
conformal metric, investigating the corresponding equations of motion in two and in four dimensions. We end the 
work in Sec. IIV1 where we include some final considerations. 

II. EQUATIONS OF MOTION AND CONSERVED CHARGES 

Let us lay down the notation, starting with the F(R) action 

S=- Jd D xVW\F(R) (1) 

where F(R) is equal to R + 2A for the usual Einstein-Hilbert action with a cosmological constant A. The equations 
of motion for this theory look like 

\ 9 acF{R) - F'(R)R ac - F"(R)[V a W c R - V b V b Rg ac ] - F"> (R)[V a RW c R - g ac W d RW d R] = (2) 

Let us now compute the conserved charge associated with this Lagrangian, following the general procedure of [I2I ] . 
From above, the boundary term one obtains from the Lagrangian is 

6 b = -F'(i?)(V Q (<5 g Qb ) - V b (6g)) - V a F'(R) Sg ab + V b F'(R)(Sg) (3) 

that is to say, the full variation of the action can be written as bS£ — \f\g\F a b Sg ab + \f\g\^aO a , where E a b are 
the equations of motion (J2J). Now, the Noether procedure gives us a conserved current associated with a symmetry 
generated by a vector field £ a . Those types of diffcomorphisms change the metric as Sg a b — V a £b + Vf,£ a - The 
conserved current is J a [£ b ] = a — £ a Jz?.. The existence of this current stems directly from the proof of the theorem. 
The hard step in any theory which is invariant under diffeomorphisms is to find the charge associated with J a . 
Consider the Poincare dual to J a : J Ql ... arl _ 1 = e a ai...a„_i^ a , where e 0l ,.. On is the volume form associated with the 
metric. We will exceptionally use n for the dimension of space-time in this subsection to avoid cluttering the notation. 
Conservation of the current implies that, on shell, the n — 1 form J is closed dJ = 0, so it can be written locally as 
a exterior derivative: J ai ... a „_i = V[ ai Q a2 ... 0rl _j = d[a 1 Qa 2 ...a n _ 1 ], or, analogously with the duals J b — V a Q' ^. We 
can check directly that this equation above implies the conservation of the current 

V b J b = V fc V a Q[ Qb l = -i[V a ,V b ]Q[ afc ] = ^(Rabc a Q [cb] + Rabc b Q [ac] ) = \{-RbcQ [cb] + RacQ [ac] ) = Rac.Q [ac] = (4) 

If we are able to find such Q ab we can write directly the conserved charge by integrating the corresponding n — 2 
form Q over a suitable n — 2 hypersurface. The procedure is then analogous to the definition of conserved charges in 
field theory which can be written under certain assumptions as the integral of a charge flux over a sphere (a n — 2 
dimensional hypersurface in K 3 ' 1 ) "at infinity". 

As it turns out, one can indeed write a Q ab for the theory above; we suppose it has the form Q ab {{, c ] — 2W^ a ^ + 
2XV^ a ^ b \ for some metric-dependent W a and X. Taking its divergence, we will find several terms 

v a Q ab = i b v a w a - cv a w b + w a v a i b - w 6 v,£° + V a x(v a ^ - v^°) + x(y a v a t b ~ v Q v b C) (5) 

On the other hand, the direct calculation from (JTTJ) gives 

J b = F'{R)[V a {V a £ b + V b C a ) - 2V b V a C a ] ~ F"(R)V a R(V a £ b + V b £ a - 2g ab V c £ c ) - 2F'(R)V b RV a C + £ b F(R) (6) 

Comparing the terms involving two derivatives of £°, we find that X must be X — F'(R), up to some term proportional 
to the equations of motion. The terms involving only one derivative of £ a 

W a \7 a £ b - W b V a C + V a AT(V a £ b - V b D = -F"{R)V a R[\7 a Z b + V b £ a - 2 3 ab V c £ c ] (7) 
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which gives Wb = —2F"(R)VbR- We now compare the term proportional to £ a in (O and ([6]) to get 

t b V a W a - £ c V c W b = -2F'(R)R b d £ d + F(R)(, d = {F(R)g ab - 2F'{R)R ab )£ a (8) 

With the choices above for X and Wb one has, for the left-hand side 

£, b V a W a - £ c V c W b = -2F"'(R)Z b V a RV a R - 2F"(R)£ b V 2 R + 2F'"(R)C^ a RV b R + 2F"(R)£ a V a V b R 

= [2F"\R){V a RV b R - g ab V c RV c R) + 2F"{R)(V a V b R - g ab V c V c R)]t a ' ' ' 

which is the left-hand side of the equation ijHJ) by the equations of motion The conserved charge is thus given by 
Qab _ 2F'(R)V^ a ^ + 4F" (R)^ a S7 b ^ R. In this expression the first term is analogous to the Komar mass term in usual 
Einstein-Hilbert whereas the second term is new. 

Now we apply the expression for the charge found above to spherically symmetric spaces. As all of the solutions 
discussed here are static, the total mass will be given by the expression above when the vector field associated 
is the time translation vector £ a = (d t ) a , which is hypersurface orthogonal and its derivative is given by V a ^ b = 
— £[aVfc] log |£ c £ c | = (dt)\n(dr)h] dr.f(r) = e a bf'(r), where in the first step we used the general expression found in 
equation (C.3.12) from [l3j, the second step uses the fact that g a b£, b = f( r )(dt) a , and the third step defines the 
volume element of the "t — r" plane as e a b — {dt)\ a (dr)a. Last but not least, we have used the chain rule to 
expand V a f(r) = —(dx) a f'(r). One can also derive this result by noting that the relevant part of the metric 
for spherically symmetric spaces is given by ds 2 _ t = f(r)(dt) 2 — (dr) 2 /f(r). Notice that the expression for the 
derivative of £ a above allows us to write the charge as a term proportional to the spatial volume element £ a e a ai---a.n-i ■ 
Qai-a n - 2 = -2(F'(i?)V b log |£ c £ c | - 2F"( J R)V b i?)<f l e Qf)ai ... Q „_ 2 . Applying the same procedure for V a i?, we can write 
the charge as 

Q ab - 2(F'(R) f(r) - 2F"{R) R')e ah . (10) 

We note that, since the space is spherically symmetric, the integral over a D — 2 dimensional sphere one has to perform 
to compute the mass is trivial. One should remark that this expression is really better suited for asymptotically flat 
space-times. The definition of mass on non-asymptotically flat spaces is murky, where one usually uses the trick of 
"shifting" the value of the mass by subtracting the "vacuum" value, i. e., the value of Q on a maximally symmetric 
solution with the same asymptotics. We would to conclude by commenting that flat-space-solutions like Schwarzschild- 
de Sitter will have its mass modified by the presence of the F'(R) term, although not by much. One does expect, 
however, that solutions with long-range modifications of the curvature, such as those with R « 1/r asymptotically, 
will show a dramatic effect, specially for F(R) involving 1/R terms, for then F"(R) can be quite large. 

Let us now focus attention on stability. We investigate small deviations of an Einsteinian space solution R a b = %9ab 
of the general equation of motion ([2]). In order to do this, let us take the variations of the Ricci tensor and the scalar 
curvature and write them in terms of the linearized metric components h a b = Sg a b- 

SR ab = \^ 2 Kc + ^V a V c /i - V b V (a h c) b = ^V 2 h ac - -^_L_v 2 / l5ac - V b V( a h c ) b (11) 

and hence 

SR = -h ab R ab + g ac 6R ac = 2R - h — V 2 h. (12) 

y ac D ( D _ 2 ) D-2 y ' 

In the expressions above, h a b = h a b — \hg a b- Such a choice stems from the gauge invariance of linearized gravity: h a b 



and h a b + V Q «fc + ^ b v a correspond to the same perturbation 14| . By writing the variations in this form one can get 



rid of the first term in (jlip by choosing v a such that: V 2 v a + R a b Vb = —'Vbh a b . This fixes some of the gauge symmetry, 
but the equation above determines v a up to a vector w a which satisfies V 2 w a + R a b Wb = 0. We should point out that 
solving the two equations above amount to find a vector potencial in Lorentz gauge for some distribution of current. 
With this choice, we can write the first term in (fTTj) as 

VfcV^)" = V (Q V| b |/i c) b + [V h , V (a ]h c) h = -Rf ca)b hd b + R b db( ahc) d - -R d cab hd + R d ( a K) d (13) 
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We can now turn to the variation of the equations of motion ((2|) . The relevant terms are 

\h ac F{R) + ± 9ac F'(R)SR - F"(R)6RR ac - F'{R)5R ac - F"{R)[V a V c 5R - V 2 6Rg ac ] = (14) 

and we omitted terms which involve derivatives of R, which are zero on shell since R is assumed constant. One 
can further simplify the equation above by noting that the variation of the curvature scalar only depends on h, or, 
equivalently h, the trace of the metric variation. By taking the trace of the equation above, on the absence of matter, 
one finds a fourth order equation for h, which can yield the trivial solution h = if one picks the right initial values 
for the lower derivatives of h. This can be accomplished by a suitable gauge transformation for w a . So, in short, we 
can, in a general F(R) theory, accomplish the transverse traceless gauge for metric perturbations: V Q /i ab = 0, h = 0. 
In this case h a b = h ab and then the only relevant terms in (|14p are those that depend on h ab . Picking those, we have 
a simple equation for the metric perturbations 



\F{R)h ac - F'(R) Qv 2 /i oc - R b ac d h bd + R d[a h c) d ) = ( I V) 



or, specializing to Einstein spaces, 

V 2 /w + 2C b ac d h bd + j^jh ac - §7^h ac = 0, (16) 

where C a bcd is the Weyl tensor. The causal structure of the equation above for Einstein spaces has been studied in 
15 1, and we will take the result of a particularly nice survey of [lq . As one can see from, for instance, the behavior 



of the equation above under conformal transformations, the causal structure of the equation above is governed by a 
"mass term" which depends on the spin of the particle as well as the scalar curvature. In this case, the effect is that 
the mass term M 2 of spin 2 fields is modified by the presence of curvature to 



T 2_ F{R) , 2 



m =-fW) + d r (17) 

Note that this means that the term proportional to the Riemann tensor in ([15]) takes into account the causal structure 
of spin-2 particles and hence only the following terms contribute to the mass of the perturbations. This is essentially the 
spin-2 version of the Breitenlohner-Freedman found [It]] in the study of anti-de Sitter spaces. One can also check that 
constant curvature solutions of the generic Lagrangian (which includes Einstein- Hilbert with a cosmological constant 
as a particular case) always have massless perturbations (gravitons). Therefore, these solutions do not present us 
with new instabilities. These are essentially the same which arise in the study of Einstein-Hilbert Lagrangian. 

The equation above (TIT)) gives a straightforward criterion for the stability of the solutions described below: if the 
effective squared mass of the gravitons is positive, small deviations of the metrics found will oscillate, otherwise they 
will grow exponentially and the space found will be unstable. For the special case D — 4, the corresponding equation 



has been found by several authors - see [18[ and references therein. 

Theories involving 1/R terms, which we call "non-perturbative," has attracted some attention recently, mainly 
because of the prospect of having a solution of positive curvature, i.e., a assymptotically de Sitter space-time, without 
the need for a cosmological constant. We analyze one well-known solution of those theories which shows an instability 
against curvature perturbations which should plague those theories in general. 

We remark that the mere existence of a term proportional to the inverse of the curvature scalar does not fit 
with what we know from renormalizable field theory, or even what little we can grasp of quantum gravity, as in 
string theories. Usual corrections from perturbation series in the interactions of quantum fields with classical gravity 



show the dynamical generation of terms up to R D I 2 [19]. In string theory, general remarks given in 20 1 can be 
applied here to constrain what corrections one can expect from string instantons such as those that arise from 
compactifications which respect N = 1 SUGRA. Such compactifications have received a lot of attention in recent 
years since the work of 21] introduced extra ingredients like confined fluxes, which allowed for meta-stable de Sitter 
vacua. At any rate, instantons and fluxes may break the Super- Weyl rescaling symmetry of SUGRA \g ab — > e 2A <? a &j 
and e 2A = — 3/(<f>^"$i + Ci&i + Ci&f — 3), where $i is a chiral superfield] to a discrete set, as it happens with other 
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superfields. The residual symmetry constrain the effective potential generated to be of the form e~ aR . In this 
expression, a incorporates scalar condensates, but should not depend on the metric. One should note that even in 
this hypotetical scenario, the corrections one expects from quantum strings are perturbative in R. 

All the preceding discussion may present an obstacle for the consideration of non-perturbative theories, and below 
we will present a concrete criticism, following the study of [7j, while keeping the dimension arbitrary. We consider 
the general equation of motion ([2]) for F(R) = R — 7 2 /-R. For R = Rq constant, the trace of the equation gives an 
algebraic condition: 

?(*-£)-K)« 

which gives a small curvature for 7 small, as desired. We will be interested on the positive curvature solutions. Now 
consider a metric perturbation such that the curvature scalar is no longer a constant. Taking the trace of @ we find 
that 

— (D — 1){F"V 2 R + F'"(VR) 2 ) = -(D - l)V 2 (F'{R)) = yF(fl) - F'(R)R. (19) 

Taking a small curvature perturbation R = Rq + 5R, we will have 

(D - l)F"(i? )V 2 6R « - f^^F'(R ) - F"{R Q )R Q \ 5R (20) 

where we omitted a term involving one derivative of SR, which depends on the particular metric perturbation but 
will not change the conclusion. The equation above gives an effective mass for such perturbations 

m ^-^-R 0+ D -\ F ' iR °\ (21) 
D-l 2(D - 1) F"(Ro) v ' 

Note that, if we have the Einstein-Hilbert term, F'(Ro) rj 1 + 0(Rq). Now, for perturbative theories, F"(Rq) rj 1 
is independent of R and then the last term of the right-hand side dominates for small curvatures. However, for 
non-perturbative theories like F(R) = R — 7 2 /i? we will have problems; indeed, in this case we get 

9 D + 2 , x 

which shows that such theories are generically unstable under curvature perturbations. We note that this result is 
below the negative mass bound since Rq > 0. One may wonder, like in Q, whether higher order terms, perturbative 
in R, may help the situation. From the discussion above, one sees that if a term like — 7 2 /i? is present, then 
F"(R) = —2j 2 /R 3 + O(R). For the (small) constant curvature solution, R oc 7, and then F"(R) is dominated by 
the non-perturbative term, while the first derivative of F(R) is still of order 7 , hence the expression (|2~T1) should not 
change considerably as higher order terms in R are added to the Lagrangian. One should also keep in mind that the 
same arguments given above against the naturalness of a non-perturbative term to the Lagrangian also point to the 
fact that the numerical coefficients in from of higher order terms are too small to have any dramatic effect in the 
dicussion above, as can be seen in [19|. 



III. THE DILATONIC SECTOR 



Let us now study the conformally flat space 0], where the dynamics is restricted to the conformal sector and the 
metric tensor g^uix) is given by g tiy {x) = e 17 ^^^, with a(x) standing for the conformal factor, the dilaton field. In 
this case, the scalar curvature has the form 



R 



D - 1 



2Ua 



D-2 



d a ad a a 



(23) 
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The standard Einstein-Hilbert action which includes the cosmological constant is given by S = — J d D xy/\g\(R + 2A), 
where we are using 8ttG = 1. The action modificated by the change R — -> F(R) has now the form |T]). More explicitly, 
we can write 



S = - / d u xe^ a F 



D - 1 



2Da + I ^— — \ d a ad a a 



(24) 



This form shows cleary that for D = 2, and F(R) = R + 2A, the action reduces to total derivative, whereas for D = 1 
the action is reduced to the cosmological term alone. 

The equation of motion looks like De^ a F(R) + 2e% a F' (R)j^ = 0, where prime represents derivative, that is, 
F'(R) = dF/dR. In more manifest form, after substitution of the expression (|23p for the scalar curvature R, we have 

De^F(R) - (D - l)F'{R)e^- 1 ^ [2Q7 + d a ad a a\ + 

+ 2{D - l)n(e^-^ a F'(R)) — (D - 1)(D - 2)d a (e^- 1 '> a F'(R)d a a) = 0. (25) 

This is the equation of motion we have to deal with to find nonperturbative solutions in dilaton gravity, modified 
to include the change R — > F(R). Due to intrinsic difficulties, in the following we first investigate the simpler D = 2 
case. The strategy used in the following will then be extended to the case D = 4, which is treated afterwards. 

The case D = 2 is motivated by the well-known fact that gravity in (1,1) dimensions has the direct advantage 
of eliminating intricate structures which are present in higher dimensions, easing the investigation due to technical 
simplifications specific of the (1, 1) dimensions. From the equation of motion for D arbitrary, we see that in the D = 2 
case the equation simplifies to give 

e a F(R) ~ F\R)Uu + DF'(R) = 0. (26) 



We follow the first-order formalism of Ref. [llj, so we concentrate on the presence of static solutions, searching for a 
first-order equation connecting the derivative of the dilaton with some specific function of it. The static dilaton only 
depends on x = x . Using the signature (+, — ) for the two-dimensional case, we can write the scalar curvature as 

R = -e-V, (27) 

where we are using the notation a' = da/dx, a" = d 2 <j/dx 2 , etc. In this case, the equation of motion ([26]) takes the 
form e a F(R) + F'(R)a" - d 2 F'(R)/dx 2 = 0. Direct calculation leads to d 2 F'(R)/dx 2 = F"'(R) ■ (aV' - cr"') 2 e" 2CT + 
F"{R) ■ e~ a [a" 2 + 2a' a"' — a' 2 a" — a"") , and so the equation of motion gets to the new form 

e a F{R) + F'(R)a" - F"'(R) ■ {a' a" - a"') 2 e r 2a - F"(R) ■ e" CT (a" 2 + 2a' a'" - a' 2 a" - a"") = 0. (28) 

We now proceed with solving this equation. Our first strategy is to choose F(R) in the form F(R) = R + jR n + 2A, 
with 7 being a constant. Our second strategy for solving the Eq. (f2"B"|) is to choose the ansatz 

a' = -2Ae aa/2 , (29) 

This is the first-order equation we need: it is inspired in [ll| and the specific form follows from the particular profile 
of the equation of motion (f2"5)) . The first-order equation allows writing a" = 2aA 2 e aa , a'" — —4a 2 A 3 e 3aa ^ 2 , a"" = 
l2a 3 A i e 2aa . Therefore, we find the following expressions: R = -2a A 2 e^-^ ', and a' a" - a'" = AA 3 a{a - l)e 3aCT / 2 , 
and a" 2 + 2a' a'" - a' 2 a" - a"" = -AaA A (a - l)(3a - 2)e 2a<T . Substituting these expressions into Eq. ([28]) leads to 
the simpler equation 

(9 i \ OA 
1 _ - n ( n - 2)(a - l) 2 - -n{a - 1) (3a - 2) = 0. (30) 
a a J 1 

which can be examined in the two distinct cases A ^ and A = 0. The cases n — and n = 1 are trivial and will not 
be considered. 
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If A ^ 0, we have to make a = 1. Thus, for a = 1 the equation ([30]) reduces to the following equation for 
the value of A: (-2A 2 ) n (n - 1) - 2A/ 7 = 0. For n odd we get A = ±\\J - {2A/(j(n - l))) 1/n , and for n even, 
A = ±i 2; \/2A/(7(n — 1)). These values of A should be substituted to the ansatz equation (|29|) which yields the 
solution 

a(x) = -2\n(Ax + C), (31) 

where C is an integration constant. This is similar to the Liouville-like kink solution found in Ref. [221 ]. We note that, 
for a = 1 the curvature given above is constant and negative, equal to R = —A 2 ; thus, we have just found an anti de 
Sitter (adS) solution. 

Note also that the condition of reality of the dilaton field cr(x) restricts possible signs of A and 7, and the situation 
can differ for different values of n. For example, in the case n = — 1 we get 

a(x) = -21n( ±l x flx + c) (32) 



. " 2 V A 

As a result, we find that both the constants A and 7 should be either positive or negative. In the other interesting 
case, n = 2, one finds that A = ± ^/A/j. 

If A = 0, we do not need to choose a = 1 anymore; instead, the equation (|30p is now reduced to the following 
equation relating a and n: 1 — (2n/a)(n — 2) (a — l) 2 — (3/a)(a — 1) (a — 2/3) = 0. The solutions of this equation are 
a± = 2n/(2n — 1), and 0,2 = (n — l)/n. The ansatz equation in this case yields 

a(x) = -~\n(Aax + C). (33) 
a 

Again, we get to similar Liouville-like kink solution. However, in this case the value of A cannot be fixed. We notice 
that the n = 1/2 case is very specific, yielding the single root a = — 1. Also, we note that this solution does not 
depend on 7. 

Returning to the cases n = 2 and n = — 1, we find that in the A = situation both these case are possible yielding 
for n = — 1 a = 2/3 and a = 2 (for instance, for a = 2 the solution looks like a(x) = — ln(2j4a; + C)) whereas for 
n = 2 one finds a\ = 4/3 and ai = 1/2. We note that in this case, the constant curvature is not a solution anymore. 
The curvature is R — —aA 2 e^ a ~ 1 ' > ' y = — aA 2 (Aax + C) °~ 2 , so one can see that for a > 1 the curvature decays to zero 
as |cc| increases to large values. We find that for n — 2 the localized solution gives a = 4/3, whereas for n = — 1 one 
obtains a = 2. 

Perhaps the most interesting result is that the value A = somehow induces a phase transition: although the 
solutions do not change appreciably, the scalar curvature cannot be constant in the limit of a vanishing cosmological 
constant. 

Before ending this section, let us examine the possibility of finding adS solutions for F(R) general. Here we keep 
considering the first-order ansatz equation (|2l))) for the dilaton, thus the curvature P7| yields 

R = -2aA 2 e {a -V a (34) 

We get to constant curvature with a = 1, leading to adS geometry with R = —2A 2 . Indeed, we have found that the 
constant curvature R is related to the dilaton field by a" — —Re' 7 . Thus, with a = 1 the equation of motion (|28p is 
reduced to the much simpler equation 

F{R) - RF'(R) = 0. (35) 

Let us consider two more possibilities for the function F(R): (i) F(R) = R + aR n + bR m + 2A, and (ii) F(R) = 
e -aR _|_ j n t ne first case, we find that the curvature is related to the constant parameters a and b via equation 
(1 — n)aR n + (1 — m)bR m + 2A = which shows that for A ^ the flat space solution R = is impossible. In the 
second case, the equation is e~ aR (\ + aR) + 2A — 0, which admits the flat solution R = 0, for 2A = — 1. 
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We now move to four dimensional space-time. We use the same methodology, but here the equation of motion (I25|) 
takes the form 



2e 2a F(R) - F'(R)e° 



2Ua + d a ad a a 



3D(e ,T F'(i?)) - 3d a (e a F'(R)d a a) = 0. 



The scalar curvature in this case is equal to 



R = 7 



2Ua + d a ad a a 



(36) 



(37) 



Let us choose the function F(R) again in the form F(R) = R + r yR n + 2A. Thus, the equation of motion takes the 
form 



2a 



7(2 - n)R n e 2a + 3n 7 (e <T D J R 



n-l 



d a R n - 1 d a e a ) + 4Ae za = 



2a 



We substitute the value of scalar curvature (|37|) to get 



2 6 



2Da + d a ad a a 
3\ n 



+ 7 (2 -n) ( - ) e < - 2 - n ^{2D f j + d a ad a a) n + 



(38) 



(39) 



□ (e~ ( "- 1 ) <T (2Da + d a ad a a) n -^ + d a [e'^-^ (2Ua + d a ad a a) n ' 1 ^j d a e° 



4Ae : 



2a 



0. 



Here we have used □e' 7 = e' 7 (Da + d a ad a a). We recall that recently, in Ref. [23j the authors investigate cosmological 
solutions for F{R) gravity. Here, however, we search for static, spherically symmetric solutions. Thus, we need to 
replace □ by -A ~ ~d 2 /dr 2 - (2/r)d/dr, and d a Ad a B by -(dA/dr)(dB/dr) = -A'B 1 . As a result, we can write 



2Da + d a ad a a ~ -2a 1 ' 



-a - (a) 

r 



T. 



(40) 



And now the equation takes the form 
'3\ n 



a T + ^0 7 (2 - n)e {2 - n ^T n + 



2 1 1 



n 7 



dr 2 



— ( e -(™- 1 )' J T™ _1 ") — -e CT 
cir V / dr 



AAe 2a = 0. 



To solve this equation, we suppose the geometry to be adS, taking R = — (3/2)Co = const. This leads us to 



T = 2a" 



,'\2 



C e° 



(41) 



(42) 



We substitute the above equation into the equation of motion (|4"Tj) to get to the simpler equation relating Co and A: 
-3C + 27(2 - n) (-3C /2)™ + 8A = 0. The dilaton field being the solution of Eq. |42]) looks like a(r) = -2 In r, which 
makes Co to be zero. In this case, both R and A vanish. In other words, an exact spherically symmetric constant 
curvature solution is possible for A = and corresponds to the flat space. We can go further and try another solution 
of (j42|) . for Co very small. We change Co — > e, and consider a(r) — —2 ln(r) + eg(r). We get a{r) = — 2[1 + (e/4)] ln(r), 
which makes the scalar curvature very small, but now the cosmological constant does not vanish, being A = (3/8)e. 
We notice that for e positive, the solution leads to AdS geometry with A positive, and for e negative, we get to dS 
geometry with A negative. 

Another possibility in D = 4 is to use the domain wall ansatz in which the dilaton depends only on a single spatial 
coordinate, for instance a — a(x). For a review on domain walls in supergravity, see Ref. 24|. In this case we replace 
□ by —d 2 /dx 2 , 2Ua + d a ad a a by — 2a" — (a') 2 etc. The equation of motion becomes 



e°{-2a" - {a 1 ) 2 ) + 7 (2 - n)e {2 -^° (-2a" - {a') 2 ) n + 

- e°^{e-^-^(-2a" - (a') 2 )"" 1 ) - 4~ ( e^ 1 ^ {-2a" - {a 1 ) 2 )^ 1 ) e°a' 
dx l dx \ I 



(43) 



2(1 



+ 4Ae 2 °" = 0. 



9 



where prime now stands for derivative with respect to x. Here we used the fact that -^e a — e a a' . To solve this 
equation, as before we consider the ansatz a' = —2Be°l 2 . It allows writing 

2a" + (cr') 2 = 85V. (44) 

This leads to (-12S 2 )"7(2 - n) - 12B 2 + 4A = 0, and for A / we can obtain the constant B. 
We solve the equation (|44|) to get 

a(x) = -2hx(Bx + C), (45) 

with C being an integration constant. The scalar curvature in this case is R = — 12B 2 . 

If we impose the condition that the scalar curvature (|37|) is constant, we arrive at the following reduction of the 
Eq. (|36|) : 2F(R) — RF'(R) = 0. Thus, if we follow the steps done in two space-time dimensions, we can consider other 
possibilities for the function F(R): (i) F(R) = R + aR n + bR m + 2A, and (ii) F(R) = e - aR + 2A. In the first case, the 
curvature turns out to be related to the constant parameters a and b via equation R+(2 — n)aR n + (2 — m)bR m +AA = 
which shows that for A ^ the flat space solution R = is impossible. In the second case, the equation is 
e~ aR (2 + aR) + 4A = 0, which leads to vanishing curvature for 2A = —1. Note that, by the criterion (fT7|) . such 
solutions are always stable for n, m > 0. 



IV. ENDING COMMENTS 

In this work we have studied general features of F(R) gravity models and derived a formula for the asymptotic mass 
of such theories, showing that for perturbative theories in which F(R) allow a perturbative expansion the value of 
the mass do not change significantly from that of usual Einstein-Hilbert solutions. We have analyzed the behavior of 
these theories under perturbations of the solutions and concluded that, while metric perturbations are essentially the 
same as in Einstein-Hilbert models, the scalar curvature perturbations are unstable for positive curvature solutions 
of models which modify the usual Lagrangian by R n terms with n < 0. It would seem that, on top of providing 
an effective action mechanism by which such term would arise, one is riddled by the stability of the solution found. 
Static solutions of the equations of motion are found for the F(R) gravity models restricted to the conformal sector. 
We see that these solutions have the Liouville-like kink form in two space-time dimensions. Also, very similar effects 
are achieved in four-dimensional space-time, both for spherically symmetric or domain-wall-like solutions. We also 
have found that the limit A — > takes these solutions with R ^ 0, then still maintaining the desired behaviour of 
F(R) theories to present curved backgrounds without the need for a cosmological constant. The presence of first-order 
equation which solves the equation of motion with adS geometry indicates that the system may admit supersymmetric 
extensions, at least for n positive. If this is indeed the case, the supersymmetric realization will be new, given the 
non-standard form of the graviton Lagrangian. This is an issue in which we will report elsewhere. 
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